Imaging Studies of Characteristics for a Slab of a Lossy Left-handed Material 
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The characteristics of an imaging system formed by a slab of a lossy left-handed material (LHM) 
are studied. The transfer function of the LHM imaging system is written in an appropriate product 
form with each term having a clear physical interpretation. A tiny loss of the LHM may suppress 
the transmission of evanescent waves through the LHM slab and this is explained physically. An 
analytical expression for the resolution of the imaging system is derived. It is shown that it is 
impossible to make a subwavelength imaging by using a realistic LHM imaging system unless the 
LHM slab is much thinner than the wavelength. 
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I. INTRODUCTION 



Materials that have simultaneously a negative electric 
permittivity and a negative magneticpermeability have 
attracted much attention recently [E S S 13 ■ These ma- 
terials have extraordinary properties, such as the neg- 
ative refraction, inverse Doppler shift and backward- 
directed Cherenkov radiation cone, etc 0- The electric 
field, the magnetic field and the wave vector of an elec- 
tromagnetic wave propagating in such a material obey 
the left-hand rule and thus it is called the left-handed 
material (LHM). Due to its effect of the negative refrac- 
tion, a LHM planar slab can act as a lens, and focusing 
waves from a point source i . Recently, Pendry predicted 
that a LHM slab can amplify evanescent waves and thus 
behaves like a perfect lens [(j. It is well-known that a 
conventional imaging system is limited in resolution by 
the wavelength of the used light, since the information 
finer than the wavelength is carried by evanescent waves, 
which are drastically decayed before they reach the im- 
age plane. A LHM slab appears attractive in provideing 
an effective approach to achieve a subwavelength imag- 
ing. However, the LHM that Pendryconsidered is quite 
ideal and has no loss. Garcia et al 0] pointed out later 
that an actual LHM should be lossy and the LHM loss 
may completely suppress the amplification of evanescent 
waves, and consequently a LHM slab can not act as a 
perfect lens. It is interesting to find out answers to the 
following questions: (1) What is the characteristic of a 
realistic LHM imaging system? (2) How does the reso- 
lution of the system depend on the LHM loss and other 
slab parameters? (3) Is it possible to use a realistic LHM 
imaging system to make a subwavelength imaging? If 
yes, then how to realize it? In the present paper, we will 
study these problems theoretically. 



II. TRANSFER FUNCTION OF A LHM 
IMAGING SYSTEM 

Consider a LHM slab imaging system with co-ordinates 
as shown in Fig. 1. We assume that the system is 
invariant along the y-axis direction and only the one- 
dimensional imaging will be considered here for simplic- 
ity. In the imaging system, the space can be divided into 
three regions: region 1 (z < 0), where the objective plane 
is located; region 2 (the LHM slab); region 3 (z > d), 
where the image plane is located. The electric permit- 
tivity and magnetic permeability in the three regions are 
denoted by e/ and /U; (I — 1,2,3), respectively. Here 
e\ = £3 > 0, /ii = (X3 > 0, Ree 2 < 0, and Refi 2 < 0. This 
imaging system has the properties of being linear and in- 
variant under translation along the x direction, and thus 
can be described by a transfer function (Ti.) or impulse 
response (h) (the latter is just the Fourier transform of 
the former). 

In order to calculate the transfer function of the LHM 
imaging system, we consider a plane wave propagating 
along the z direction from the objective plane z = —a 
(a > 0). The plane wave will refract and reflect simulta- 
neously at the interfaces z = and z — d before it reaches 
the image plane z — d + b (at a distance of b away from 
the LHM slab). The electromagnetic fields can be de- 
composed into the E- and 77-polarization modes. Let's 
first consider the _E-polarization case, where the electric 
field has only y component (i.e., E — Ey), and the mag- 
netic field has the form of H = H x x + H z z. All the field 
components are assumed to have the time dependence of 
form exp(— iuit). The electric field E of the wave in the 
three regions can be written in the following form of 

E\(x,z) = exp(ik x x) {Ai exp(ik z (z + a)) + B\ exp(— ik z zj^) 
E 2 (x,z) = exp(ik x x){A 2 exp(-iqz) + B 2 exp(zgz)}, (2) 
E 3 (x,z) = A 3 exp(ik x x) ■ exp(ifc z (z - d)), (3) 
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where k x is the x-component of the wave vector, k z 
is the z-component in regions 1 and 3, and q is the 
z-component of the wave vecto r in region 2. Obvi- 
ously, one has k z — n\k^ — k% for propagating waves 
(i.e., k x < niko, where n\ — ^/ei/ii), and k z — 
i\/k^ — n\k,Q for evanescent waves (i.e., k x > n\ko). 

fc 2 for propagating waves and 
The coeffi- 



n\kl 
Similarly, q = 



for 



62^2^0 ^ 0r evanescent waves. 



q = i^fkl 

cients for the forward and backward (concerning the en- 
ergy flow) propagating components are denoted by A and 
B, respectively. Note that the energy flow for the term 
A2 exp(-iqz) in Eq.(1.2) is along the +z direction since 
the refraction index of the LHM slab is negative. 

The transfer function of the imaging system is defined 
as H.(k x ) = E 3 (x, d + b)/Ai exp(ik x x), i.e., TL(k x ) = 
(A3/A1) exp(ik z b). From Maxwell's equations, one has 
H x = (i/LU[j,Qii)(dE/dz). From the continuity condition 
of the field components of E and H x at the interfaces 
z = and z = d, one can derive formulas for A3/A1 and 
TL{k x ). In the present paper we can write the transfer 
function Ti.(k x ) in the following product form (with each 
term having a clear physical interpretation) 



H(k x ) 



hPt 2 P, 



where 



U = 



1 



P 
P 



2 1- 

exp(— iqd), 
exp{ifc z (a + b)} 
1 / Mig H2k z 
2\fi 2 k z Hiq 



(4) 



(5) 

(6) 

(7) 
(8) 

(9) 



Here t\ represents the effective transmission coefficient 
for a plane wave at the interface z = (i.e., t\ — 
A2/(Aiexp(ik z a))), while t 2 is the simple transmission 
coefficient at z = d (i.e., t 2 = A3 / (A 2 exp{—iqd))) . £ 
is related to the impedance matching of the LHM with 
the surrounding medium (£ = — 1 when the impedance is 
matched). The second term (1+£)P 2 in the denominator 
of expression (3.1) for t\ is attributed to the effect of the 
reflection of the wave at the interface z — d. In fact, this 
term ((1 + (,)P 2 ) vanishes if the LHM slab is completely 
matched in the impedance with the medium in region 3 
(i.e. e 2 = —€3, [I2 = — M3) regardless what the value of ei 
or /xi is. P is the phase change factor or amplitude am- 
plification factor as the wave passes the LHM slab, and 
P is the total phase change factor or amplitude decaying 
factor as the wave goes from the objective plane to the 
left surface of the LHM slab and from the right surface 
of the LHM slab to the image plane. 

For an ideal LHM imaging system (i.e., £2 = — ej., ^2 = 
—fii and b — d — a), one has £ = — 1, t\ = 1, t% = 1, and 
7~L{k x ) — 1 for all k x , this means the imaging system is 
perfect. 



III. INFLUENCE OF THE LOSS IN THE LHM 

It is known that an actual LHM should be lossy (i.e. 
e 2 or [12 has an imaginary part) since any LHM must be 
dispersive @- Here we analyze the characteristics of the 
transfer function for an imaging system formed by a lossy 
LHM slab. Consider first the case when the loss comes 
only from the imaginary part of e 2 , i.e., 

£2 = -ei(l -iS), H2 = -^i, 

where we assume < 5 <C 1. As Ti. is an even function 
of k x , we only needs to consider the interval [0, 00) for 
k x . For the convenience of analysis, we further divide 
this interval into two regions, namely, k x G [0,ni£;o) and 
k x € [nifcoj+oo). 

1. The case when k x < niko. Th is case co rre- 
sponds to propagating waves with k z = ^nfkQ — k x be- 
ing a real number. Then one has q — y/ e 2 [i 2 kQ — k^ — 
\/k z — iufk^S. 

(i) Consider the special situation when the plane waves 
are propagating nearly along the z axis, i.e., k x <C 
\/\ — 5n\ko (or equivalently, k z 3> v^^ifco)- In this case, 
one has £ = — 1 + 0(<5 2 ), which implies that the media are 
approximately matched and t\ sa l/i 2 ~ 1. The transfer 
function can then be approximated by 



H (k x ) w exp 



niko 



^n{kl-kl 



D6 



(10) 



where D = mkod. At the point k x — (corresponding 
to the plane wave normally incident on the LHM slab), 
one has H(0) = exp(— nikodS/2). This indicates that the 
small loss of the LHM makes H (0) a bit below 1. 

(ii) Then consider the other situation when k x — » niko 
(i.e., k z 0). In this situation, one finds that k z <C 
VSriiko and £ = —\/—iSniko/ (2k z ) differs from —1 com- 
pletely. Therefore, a small loss of the LHM causes a 
serious mismatch problem for these waves. The transfer 
function can be approximated by 



H(k x ) 



2y/n(k* - kl 



exp 



n\k§D& 

Obviously, TL{k x ) when k. 



-D 



)}. I") 



niko. Therefore, the 
existence of a small loss in the LHM changes the value 
of the Ti function completely (dropping to from 1; see 
Fig. 2(c) below) in the neighborhood of k x — niko. 

2. The case when k x > nikp. This case corresponds to 
evanescent waves with k z = i\/k 2 x — n\k^. Then one has 

q = W k l~ £ 2^2fco = W\ k z\ 2 + iSnjk^. 

(i) Consider first the situation when k x — » niko (i.e., 

k z — > 0). In this situation, one finds that \k z \ <gl VSniko 
and £ w — ^/i6niko/(2\k z \) differs very much from — 1. 
Thus, one sees that a small loss of the LHM causes a se- 
rious mismatch problem for these evanescent waves. One 
can then approximate the transfer function as 



H{k x ) 



2VW 



nikn 



nikoDS 



exp{rjy|(l + l )}. (12) 
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Obviously, Ti — > when k x — > nifeo- We find that the 
real part (7i r ) of 7i is far smaller than the imaginary 
part (Hi) in this case. Tii decreases as k x decreases to 
nifco- On the other hand, Tii should approach zero when 
k x niko 0. Thus, there exists a peak of Tii on the 
right-hand side of the point k x — niko (see Fig. 2(d) 
below). 

(ii) Then consider the other situation when k x 3> 
\J\ + Sniko (or equivalently, \k z \ S> ySriiko). These 
evanescent waves are very important for a subwavelength 
imaging. For this case our analysis shows that £ = 
— 1 + (niko/k z ) 4 S 2 /8, which indicates that S causes only 
a second-order perturbation to £ (the LHM slab seems 
still well-matched with the surrounding media for these 
evanescent waves). However, |P| w exp(y^k x — n^k^d) 
increases exponentially as k x increases. Thus, unlike the 
case of propagating waves, the second term (1 + £)P 2 in 
the denominator of the expression (3.1) for t\ can not 
be neglected. In other words, a tiny S may cause a sig- 
nificant change of the effective transmission coefficient 
t\ for the evanescent wave at the interface z = (note 
that 5 causes only a small perturbation to ti and P), 
and consequently leads to a quick decrease of \Ti\. For 
a slab of usual medium (with or without loss), the ef- 
fective transmission coefficient t\ must be close to 1 if 
the impedance is almost matched (regardless the value of 
k x ). However, for a LHM slab t\ can be very small (con- 
sequently the total transmission coefficient of the LHM 
slab can be very small) even the the impedance is almost 
matched. The main reason is that the evanescent wave 
"reflected" at the back interface of the LHM slab may 
be effectively amplified (before it comes out of the slab) 
and thus significantly suppress the transmission of the 
evanescent wave according to Eq. (3.1). Therefore, the 
conventional concept of impedance match is insufficient 
to describe the medium match for evanescent waves for 
a LHM slab. We introduce a critical value (k x ) for k x , 
at which ReTi « 1/2. Apparently, k x is closely related 
to the resolution of the imaging system. After a pertur- 
bation analysis for Eqs. (2)-(3.5), one can conclude that 
the critical value k x can be approximately determined 
from the condition i?e[(£ + 1)P 2 ] = 2 (physically this 
condition can be understood as a balance condition for 
the two terms in the denominator of expression (3.1) for 
ti). When k x 3> k x , the term (£ + 1)P 2 becomes dom- 
inant in the denominator of expression (3.1) for t\, and 
t\ decreases rapidly to (consequently Ti — ► 0) as k x in- 
creases. The transfer function can then be approximated 
as 



k x = nik (a bit larger than mho) as discussed before. 
Nevertheless, Tii can be approximated to zero everywhere 
(except in a very narrow region around niko] this non- 
zero region is so narrow that it can be neglected in an 
imaging system for the image in the physical space). Our 
numerical analysis shows that k x ~ 2.3nifco for this ex- 
ample. Fig. 2(a) shows that Ti decreases to zero rapidly 
as k x increases in the region k x > k x . From these figures 
one sees that the characteristics of the transfer function 
Ti (calculated from Eq. (2)) agree very well with our 
earlier analysis. 

The value of k c x is very important for an imaging sys- 
tem since k x directly relates the smallest scale in which 
information can be correctly transferred by the system. 
To see clearly what influences this parameter, we wish to 
derive an analytic formula for k x . From Pe[(£ + 1)P 2 ] = 
2, one obtains 7e~ 7 = p (i.e., In 7 — 7 = hip), where 
7 = ZV0£) 2 -n 2 fc 2 /(2nifc ) and p = DVS/4. We as- 
sume that loss parameter 6 is so small that p< 1. Then 
one has e 7 ^> 7 > 1, i.e., 7 ^> In 7 > 0. Thus it follows 
that 7 = — hip + In 7 « — In p. This approximation can 
be improved as 7 ~ — lnp + lnln(l/p). One can finally 
obtains the following analytic formula 



D 



D \dVsJ 



-1 2 



Note that D = nikod. 

Eq. (8) shows that for a given loss parameter^, k x 
decreases as the thickness of the LHM slab increases (ac- 
tually k%, is almost inverse proportional to d as one can 
see from 7 « — lnp). On the other hand, k x is almost lin- 



4/ (DVtj 



Ti(k x ) « 16 



n\kl 



n\k 2 5 



exp(— 2y fc 2 - n\ k^d), k x > 



Figs. 2(a) and 2(b) show the real and imaginary parts 
of TC(k x ) (calculated with Eq. (2)), respectively, where 
we have chosen d = A (A = 2ir/niko) and 5 = 1CP 4 . Figs. 
2(c) and 2(d) give the corresponding enlarged views of 
Ti r {k x ) and Tii(k x ) in the neighborhood of niko. From 
Figs. 2(b) and 2(d) one sees that Tii has a peak around 



early proportional to In S (since (2 /D) In 

— In (S) /D + 21n(16/P) /D and the first term is dom- 
inant over the second term). Thus, k c x does not in- 
crease noticeably even if 5 decrease by one order of 
magnitude. To realize the subwavelength imaging, e.g., 
k x > lOfcoj h is required that S < 10~ 52 when d = 2 A, 
and <5 < 2 x 10~ 25 when d = A. It is doubtful whether 
the loss of a realistic LHM could be reduced to such a 
level. Therefore, the thickness of the LHM slab should 
be smaller than the wavelength in order to realize the 
subwavelength imaging for such a system. 

From Eq. (2)-(3.5) (or Eq. (8)) one can find that 
the characteristics of the transfer function (as well as 
k x /(niko)) depend only on 5 and d/X. Fig. 3(a) gives 
the contour plot of fc°/(mfco) (calculated with Eq. (2)) as 
a function of 5 and d/X. From this figure one sees that 
k c x is very sensitive to the thickness d. Our numerical 
Calculation shows that the requirement for k x = lOnifco 
(for a subwavelength imaging) is 5 = 2 x 10 -25 when 
d = A, 5 = 9 x 1CT 11 when d = A/2, and 5 = 6 x 1CT 5 
when d — A/4. Fig. 3(b) shows k x obtained from Eq.(2) 
and the analytic expression (8), respectively, as the loss 
parameter varies. One sees that they agree very well. 

Consider now the impulse response for a lossy LHM 
imaging system. The impulse response h(x) is defined 
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as the field distribution on the image plane produced 
by a point source on the objective plane. Obviously, 
it is just the Fourier transform of the transfer function 
H(k x ). Our analysis shows that H(k x ) is approximately 
a step function, i.e., H(k x ) w 1 for \k x \ < k% (except the 
abrupt change in the narrow neighborhood of k x = n\ko, 
which is narrow enough to be neglected in calculating 
the impulse response by the Fourier transform), while 
TL{k x ) w elsewhere. Thus, we have approximately 
h(x) = 2k x sin(k x x) / (k x x) . According to the definition 
of the resolution (denoted by A) for an imaging system, 
we find that A = Tr/k x . Figs. 4(a) and 4(b) show the real 
and imaginary parts of the impulse response (normalized 
by h(0)) of a LHM imaging system, respectively, with 
S = 10~ 4 and d = A. One sees that the imaginary part of 
the impulse response is very small and can be neglected 
as compared with the real part. From Fig. 4(a) one sees 
that the real part of h calculated from the Fourier trans- 
form of Eq.(2) (the solid line) agrees very well with that 
from our analytic expression (the dotted line) in the re- 
gion of | a; | < 2^jk x (note that our analytic expression 
gives a real value of h; the imaginary part of h is always 
very small, as show in Fig. 4(b)). This indicates that the 
resolution expression A = ir/k x is quite reliable. Thus, 
one can use the following explicit resolution for a lossy 
LHM imaging system 



A 



D \dVsJ 



+ 1 



(— 

\nik 



In the above analysis, we have assumed that the loss 
of the LHM comes from the imaginary part of €2- We 
can also discuss the case in which ^2 gives rise to the 
loss, i.e. 62 = — ei, ^2 = — A*i (1 — iS). Our analy- 
sis and numerical calculation have shown that the char- 
acteristics of the transfer function for this case is sim- 
ilar to the previous case that we have discussed. In 

this case one has k x w (nifco) \J hi (2/S)] 2 + 1 and 

A w (n/mko) /\J[{l/D) In {2/5)f + 1. For the same d 
and d, k x for this case is a bit smaller than that for the 
previous case (see Fig. 5). 

In the above analysis the electromagnetic wave is as- 
sumed to be E— polarized. For the H— polarized case, 
one can obtain the transfer function from Eq.(2) di- 
rectly by interchanging e and \x due to the symmetry 
of Maxwell's equations. Thus, similar analytic formulas 
for k% and A can be found. 



and evanescent waves; (2) TL decreases to zero rapidly 
when k x increases from a critical k x (at which TL = 1/2). 
Our analysis shows that k x is almost linearly propor- 
tional to the logarithm of the loss parameter 6 (related 
to the imaginary part of e 2 or fi 2 ) of the LHM. Thus, k x 
does not increase noticeably even if the loss decreases by 
one order of magnitude. Analytical formulas for k% and 
the resolution of the imaging system have been derived. 
It has been shown that the thickness of the LHM slab 
should be smaller than the wavelength in order to realize 
the subwavelength imaging for such a system. 

Our analysis has also shown that a tiny loss of the 
LHM may suppress the transmission of evanescent waves 
through the LHM slab due to the amplification of the 
evanescent waves in the LHM slab. Such a suppression 
can be effectively relaxed by reducing the thickness of 
the LHM slab. We have also tried to eliminate or reduce 
the destructive influence of the LHM loss to the imaging 
resolution through manupilating the real and imaginary 
parts of H2 and £2- However, we found this is not possible 
(see the appendix). 
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Appendix: Possibility of reducing the destructive 
influence of the LHM loss of /U2 and e-2 

In this appendix, we study the possibility of eliminat- 
ing or reducing the destructive influence of the LHM loss 
to the image through manupilating the real and imagi- 
nary parts of H2 and 62 . Consider the case when the elec- 
tric permittivity and magnetic permeability of the LHM 
slab have the following forms 

£2 = -£i(l ~S), 
M2 = - 5"), 

where S = S' r + i5 i , 6 = S r + id i (S' r , S i: S r , and are 
all real). We assume that S i ,S i > 0, and \5 |, \5 | <C 1. In 
the following analysis for the LHM slab, our attention will 
be focused on the evanescent waves for the subwavelength 
imaging. 

From Eq. (3.5), one has 



IV. CONCLUSION AND DISCUSSION 

In this paper we have analyzed the characteristics of 
the transfer function for an imaging system formed by a 
lossy LHM slab. The small loss in the LHM may com- 
pletely change the transfer function in the following as- 
pects: (1) the transfer function drops to zero abruptly 
near the turning point k x — nifcg for the propagating 



1 



where <5 = (S' +S")/2, (3 = ^Jk 2 x - npj/(nifc ) (k x > 
niko for the evanescent waves). Then the transfer func- 
tion becomes 



n 



1 



■0.25(5" +<5//3) 2 exp(2D/3) 



5 



Therefore, the requirement for reducing significantly the 
influence of the LHM loss is 

d" = 0, 5 = 0. 

Obviously the above requirement is equivalent to that of 



S = and S — 0, which in fact corresponds to the case 
of an ideal LHM without any loss. Thus, it is impossible 
to effectively reduce the destructive influence of the LHM 
loss to the subwavelength imaging by manupilating the 
real and imaginary parts of 112 and ti . 
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FIG. 1: Schematic diagram for an imaging system formed by 
a LHM slab. 
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FIG. 3: (a) Contour plot of k%/{n\ko) as a function of the 
loss parameter 8 and the thickness d. (b) The dependence of 
the critical k%/(niko) on the loss parameter 8 when d = A 
or 0.5A. The solid lines and the marks correspond to the 
results obtained from Eq. (2) and our analytic formula (8), 
respectively. 
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FIG. 4: (a)The real part of the normalized impulse response 
calculated from the Fourier transform of Eq. (2) (solid line) 
and our approximate analytic formula (dotted line). (b)The 
imaginary part of the normalized impulse response calculated 
from the Fourier transform of Eq. (2). Here d — A and 
<5=1CT 4 . 
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FIG. 5: Comparison of the critical values of k% when the 
loss comes from the imaginary part of e (solid lines) and the 
imaginary part of p (dashed lines). 



